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Abstract 

The functional-integral quantization of non-Abelian gauge theories is affected by the Gribov 
problem at non-perturbative level: the requirement of preserving the supplementary conditions 
under gauge transformations leads to a non-linear differential equation, and the various solutions 
of such a non-linear equation represent different gauge configurations known as Gribov copies. 
Their occurrence (lack of global cross-sections from the point of view of differential geometry) 
is called Gribov ambiguity, and is here presented within the framework of a global approach to 
quantum field theory. We first give a simple (standard) example for the SU(2) group and spher- 
ically symmetric potentials, then we discuss this phenomenon in general relativity, and recent 
developments, including lattice calculations. 
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I. INTRODUCTION 

In the modern geometrical formulation of field theories one may start with general prop- 
erties of classical dynamical systems described globally by a trajectory or history [1]. A 
history is a section of a fibre bundle E having the space-time manifold (M, g) as its base 
space. The typical fibre is known as configuration space and can be denoted by C. If the 
dynamical system consists of a set of fields, such a set takes its values in C. One usually 
assumes that a fibre bundle E is associated with a group Q; such a group can act globally on 
the fields, or instead it can be a Lie group that acts locally on each fibre. In the latter case, 
it is more appropriate to consider a principal fibre bundle P [2] which has Q itself as the 
typical fibre, and a connection on P. The connection on P defines a connection 1-form on 
the space-time (M,g), known as a gauge field and taking values in the Lie algebra L(Q) of 
Q. The connection is therefore a collection of sections of trivial bundles C/j x T*(M) x L(Q), 
which are related to each other by group actions in the overlaps Ui fl Uj. These group 
actions are called gauge transformations, and for the whole dynamics to remain invariant 
under gauge transformations the gauge field itself has to become dynamical. The space 
<3> of all possible field histories includes both those that do and those that do not satisfy 
the Euler-Lagrange equations (when there is an action principle for the former), and can 
be viewed as an infinite-dimensional manifold. To sum up, the space of histories $ is a 
principal fibre bundle having the infinite-dimensional Lie group Q, the gauge group, as its 
typical fibre, while physics takes place on the base space of this bundle, the latter being the 
quotient space $/C? called the space of orbits [1]. 

On the other hand, the functional-integral quantization of gauge theories of fundamental 
interactions relies upon the following properties. For a given gauge theory with action func- 
tional S and generators of infinitesimal gauge transformations R l a with associated linearly 
independent vector fields 1Z a = FL l a j-j, one has (hereafter, Greek indices from the beginning 
of the alphabet are Lie-algebra indices, Greek indices from the middle of the alphabet are 
space-time indices, while lower case Latin indices are used for fields at a space-time point 



If all flows that leave the value of the action invariant can be expressed, at each point of $, 
as linear combinations of the 1Z a and skew fields at that point, the most general solution of 



[1]) 




(1) 
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Eq. (1) reads as 



Ha, Hp - c\ p + j—t 1 ^. 



(2) 



Hereafter, we assume that the tensor fields T a p vanish, and that the CTg are constant, in 
that they are independent of the fields (but they might depend on space-time coordinates). 
The are therefore structure constants of an infinite-dimensional Lie group. Since we 
are studying a gauge theory, what is crucial is the existence of equivalence classes under the 
action of gauge transformations. The information on such equivalence classes is encoded into 
suitable 'coordinates', say (I A ,x a ), obtained as follows. The I A are non-local junctionals 
which pick out the orbit where the field tp lies, whereas the x° L { i f) P^k ou t the particular 
point on the orbit corresponding to the field tp. They are precisely the gauge functionals 
considered by Faddeev and Popov [3]. The physical (out | in) amplitude can be expressed, 
formally, by functional integration over equivalence classes (under gauge transformations) 
of field configurations: 



It is now possible, at least formally, to re-express this abstract functional-integral formula 
in terms of the original field variables, under the assumption (not always verified) that 
the x° coordinates are globally defined. For this purpose, the integration over equivalence 
classes is made explicit by introducing the ^-integration with the help of a 5 distribution, 
say 6(x(<p) — ()■ After obtaining the Jacobian J((p) of the coordinate transformation from 
(I A ,X a ) to <f l , its functional logarithmic derivative shows that J takes the form J(<p) = 
N((p)detJ-, where a careful use of dimensional regularization can be applied to reduce N(ip) 
to a factor depending only on the non-local functionals I A [1]. Its effect is then absorbed 
into the measure over the field configurations, so that the (out | in) amplitude reads as 



At this stage one performs a Gaussian average over all gauge functionals, and denoting by p a/3 
a constant invertible matrix, and by a a and ip 13 two real- valued and independent fermionic 
fields [1], the functional integral for the (out | in) amplitude is eventually re-expressed in the 
form 



The physical predictions of the theory are x~ and p-independent. The operator JF^ is called 
the ghost operator, and the fields a a and ij>P are the corresponding ghost fields [1]. The 




(3) 




(4) 
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physical idea is due to Feynman [4], while the development of the corresponding formalism 
in quantum gravity was obtained by DeWitt [5]. However, to have a well defined functional- 
integral representation, it is more convenient to start from an Euclidean formulation. In 
particular, in the one-loop quantum theory, one considers infinitesimal gauge transforma- 
tions, and if the theory with action S is bosonic, one tries to choose x a m such a way that 
both the ghost operator and the operator on perturbations of the gauge field are of Laplace 
type [6]. The addition of the gauge- averaging (or gauge- fixing) term x a Pa/3X fB plays a crucial 
role in ensuring that both the ghost operator and the gauge-field operator have a well defined 
Green function. Such a term corresponds to the supplementary (or gauge-fixing) condition 
already occurring in the classical theory [5]. 

II. ASYMPTOTIC CONDITIONS FOR PURE YANG MILLS THEORIES 

In geometrical language, the Yang-Mills potentials A a (x) are components of the pull-back 
to space-time of the Lie-algebra-valued connection 1-form A a ^dx^ on the principal Yang- 
Mills bundle, and one has basis matrices G a for an /-dimensional irreducible representation 
of the Yang-Mills Lie algebra with structure constants fp such that 

[G a , Gp] = G-y p a p. 

The gauge potentials A a ^ and basis matrices G a can be used to define the matrices A^ = 
G a A°^, and one can write finite gauge transformations in the form 

D A ll = -D 7li D- 1 + DA li D- 1 . (6) 

This equation describes the gauge orbit for the configuration A^, i.e. the set of all gauge- 
equivalent configurations corresponding to A^. One can further say that finite little gauge 
transformations [1] are the proper subset of transformations described by Eq. (6) such that 
the matrix D can be written as 

D = e G "t a , (7) 

where the £ a are finite functions on space-time that vanish at spatial infinity. This condition 
implies that D tends to the identity therein. The matrices D are the representation matrices 
for the group Q generated by the G a . 
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Big gauge transformations are defined by Eqs. (6) and (7), supplemented by the follow- 
ing asymptotic conditions on the D matrices (the £ a functions no longer vanish at spatial 
infinity) : 

D{x) — > ±1; as x — > oo, (8) 

1; being the /-dimensional unit matrix. In this case, both the identity and minus the identity 
are necessary, because the transformation law (6) does not change on replacing D by —D. 

More precisely, the infinite-dimensional Lie group Q with structure constants C\q as in 
Eq. (2), with representation matrices expressed by Eq. (7), is the proper gauge group (its 
elements are the little gauge transformations defined above). The full gauge group is obtained 
by adjoining to the proper gauge group all transformations of the space $ of histories into 
itself, independent of field variables, which leave the action functional S un-affected and 
do not result from global symmetries. When big gauge transformations exist, both the full 
gauge group and the space of histories have disconnected components. 

Note that both big and little gauge transformations do not affect the dynamical variables 
at infinity. This leads to an investigation of the boundary conditions that should be imposed 
upon the dynamical variables. In general, the Yang-Mills fibre bundle might be non-trivial 
in spacelike directions. The possible physical implications of such a non-trivial nature have 
not yet been studied, and we shall assume hereafter that the space-time manifold is asymp- 
totically Minkowskian, and that non-trivial aspects can only occur in the time direction 

[!]■ 

On focusing on a four-dimensional space-time, one requires that the components of the 
Yang-Mills potential should have the following asymptotic behaviour [1]: 

Aa o ~ ^> A " ~ ^ (9) 

where e is an arbitrarily small positive number. Moreover, the time derivatives of Af are 
taken to satisfy the fall-off condition 

which implies that the field strength defined as [1] 

F°liv = A a v,n ~~ A °fi,u + f%y A\ (11) 
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satisfies the fall-off conditions 

F% ~ — - — F a - ~ 1 (12) 
Three properties are therefore found to hold [1]: 

(i) All physically admissible solutions of the field equations are included in the space $ of 
field histories, in particular those having non-vanishing total charge. For these solutions, 
the charge integral converges. 

(ii) The spatial integral of the energy density T 00 is finite, whether or not the field strength 
F£ u in the expression of the energy-momentum tensor T^ u satisfies the field equations. 

(iii) The spatial integral of the Yang-Mills Lagrangian is finite. 

The above conditions are compatible with Eq. (8), but add a more accurate requirement: 
one should require the conditions 

*>=Fl|~^, Ao~^. (13) 

The space <3> consists therefore of all bundle structures of field histories A°(x) satisfying the 
asymptotic conditions (9) and (10) in every patch of the fibre bundle. The gauge group G, 
which is the typical fibre of $, can be identified with the set of all matrix functions, having 
the exponential form in Eq. (7) and satisfying the asymptotic conditions (8) and (13). The 
transition functions between patches of the bundle are matrix functions of the same form, 
restricted to overlapping regions. Last, note that the vanishing trace of the generators G a 
implies that det D(x) = 1 for all x. 

In the Hamiltonian formalism, a frequently used supplementary condition is the axial 
gauge. The resulting ghost operator reads as [1] 

T%, = -S^^x') - f« p A\5(x,x'), (14) 

the Green function of which, with indices suppressed for simplicity, is 

G(x,y) = U(x\y 1 ) + (a-l)%\x 1 )l^y>(x 2 , 2 / 2 ) ( 5(a: 3 , 2 / 3 ) 



x V exp 



f UA^x^z^x^x^dz 1 
Jy 1 



(15) 



6 



where a is any real number, and V denotes path ordering [1]. Moreover, the f a are matrices 
built out of the structure constants f& , such that 

[fa , fp\ — f~i /q/J) 

and they generate the adjoint representation of the gauge group [1]. 

One of the apparently attractive features of these Green functions is their local nature in 
the time variable, i.e. ghost fields do not propagate. However, the axial gauge is incompatible 
with the asymptotic conditions (9), (10), (12), (13) which are necessary to define the space 
of field histories and the gauge group with a certain control of the limiting behaviour at 
spatial infinity. 

To prove the latter statement, let us consider the standard gauge-fixing condition x a — 0. 
One can assume that, to fulfill such a condition, one has only to perform a finite gauge 
transformation (see Eq. (6)) such that 

= D A 1 = -D A D' 1 + D A 1 IT 1 , 

or, equivalently, D t i — D A\. This equation can be solved on each patch of the fibre bundle, 
by simply integrating with the condition that D = 1; at x 1 = — oo. In general, however, D 
does not wind up equalling the identity at x 1 = oo for all a; , x 2 and x 3 , even if the potentials 
A a ^ satisfy the asymptotic conditions (9). The same occurs when the x a are se t equal to 
any other functions of the x^ variables that vanish at spatial infinity [1]. 

For a supplementary (or gauge) condition to be compatible with the asymptotic conditions 
(9), (10), (12), (13) the resulting Green functions of ghost fields should be non-local in the 
time variable. For example, the temporal gauge A a = is instead compatible with the 
asymptotic conditions given above [1]. 

III. GRIBOV PHENOMENON 

The careful assignment of asymptotic conditions on gauge potentials is justified by the 
analysis of global properties of the latter. In the case of Yang-Mills theories, a new phe- 
nomenon is found to occur with respect to quantum electrodynamics. 

Since the discovery of instantons [7], the interest in global properties led to the attempt 
of understanding quark confinement as well by virtue of such global properties [8,9], while 
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perturbative methods failed to provide any clue on such aspects of gauge theories. Among 
the peculiar features of global nature, the Gribov ambiguity deserves special attention. In 
the simplest possible terms, on requiring that the supplementary (or gauge-fixing) condition 
should be preserved under gauge transformations of the potential, one obtains a non-linear 
differential equation, for each choice of gauge-fixing compatible with the asymptotic condi- 
tions (9), (10) and (13). The key remark is now that, on the space of orbits O = which 
is the base space of <3>, the choice of gauge-fixing on the potentials: 

X[A] - C = 

describes a certain surface which intersects [10] the orbits of Eq. (6). In the very definition 
of gauge fixing at perturbative level, one requires that each orbit should intersect the gauge- 
fixing surface only once, so as to ensure uniqueness of the potentials satisfying dynamical 
equations for a given choice of gauge fixing. In the attempt of extending such a procedure 
globally all over O, however, one faces difficulties resulting from the non-trivial topology of 
the fibre bundle in the non-Abelian case [11]. Thus, the gauge-fixing surface intersects some 
orbits more than once. In other words, in the overlapping regions among various patches 
of the bundle, transition functions are not single-valued; they might remain single-valued, 
however, within each topological sector, //this is the case, the perturbative evaluation of 
in-out amplitudes through the gauge- fixed functional integral is justified. By contrast, at 
global level, transition functions may take different values in the various topological sectors, 
by virtue of particular gauge transformations which relate them. By virtue of these big 
gauge transformations, two potentials close to each other in a given patch with respect to 
the metric on the space of histories [1] 

1$? = K YmV-9 7a/3 g^S(x,x'), 

where K YM is the Yang-Mills coupling constant, 7 Q/ 3 = —tr(f a fp) is the Cartan-Killing 
metric and g^ v is the contravariant form of the space-time metric tensor, can turn out to be 
quite far apart in another patch, since one of the two can now lie in a different topological 
sector. 

Indeed, the Gribov phenomenon can be viewed as a topological obstruction to achieving 
continuity in the choice of values ( a taken by the gauge-fixing functional x a - m geometrical 
language, one can say that it is impossible to find a global cross-section for the Yang-Mills 
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fibre bundle. Singer [12] has indeed proved that, for SU (N) groups with base manifold S* 3 or 
S* 4 , no gauge-fixing condition exists which is compatible with the asymptotic conditions (9), 
(10), (12), (13). Since the more involved symmetry groups admit SU{2) as a sub-group, the 
above statement holds for all physically relevant gauge theories. It should be pointed out 
that, in the case of quantum chromodynamics, two particular gauge-fixing conditions exist 
for which the Gribov ambiguity does not seem to occur: the axial and temporal gauges. 
However, the former does not satisfy the asymptotic conditions (9), (10), (12), (13) as we 
have seen at the end of Sec. II, while in the latter the Gauss law does not automatically 
hold. If the Gauss law is also enforced, one again finds the occurrence of Gribov copies [13]. 

Along the lines of the original Gribov argument [14] , let us consider a non-Abelian theory 
with potentials subject to the infinitesimal gauge transformations 

SAl = 8^ + r iP K X 13 = X 13 . (16) 

Moreover, let us impose a covariant gauge condition such as the Lorenz gauge for Yang-Mills 
theory: 

= 0, (17) 

for potentials satisfying the asymptotic conditions (9) and (10). The general field variables 
(p l of Sec. I are now the Yang-Mills potentials A a (x), and Eq. (17) is an example of gauge 
choice implemented, in the path integral (4), through the Dirac distribution 5(x( t f) ~ 0- 
The requirement of preserving the gauge- fixing condition (17) under gauge transformations 
yields 

= 9" [A a ^ + 5A a ^ - d^A^ = d»D a ^ X 13 = T% X 13 . (18) 

This can be seen as an equation defining the singularities of the ghost operator, i.e. its 
zero-modes. More precisely, one can interpret Eq. (18) as an eigenvalue equation for the 
ghost operator with vanishing eigenvalue. If the corresponding eigenfunction is non-trivial, 
it is said to be a zero-mode, and the ghost operator becomes singular. A physics-oriented 
interpretation is instead as follows: there exist a number of gauge potentials, related to 
each other by a gauge transformation, and satisfying the dynamical equations with a given 
choice of supplementary (or gauge-fixing) condition, when the ghost fields if) 13 have zero- 
mass bound states. The latter interpretation is suggested by the fact that JF^ occurs in the 
action functional of the theory only through the term ia a J-' c Li/j 13 in Eq. (5). Equation (18) 
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is therefore equivalent to having 



T°a\A^ = e i/; a , with e = 0, 



(19) 



where the parameter X 13 has been replaced by the ghost field ip^. 

On assuming, for simplicity, that the gauge potential vanishes before performing the 
gauge transformation, Eq. (19) becomes, for a generic value e of the ghost eigenvalue, 



This equation can only be solved for non-negative values of e. On taking increasingly large 
magnitudes of the potentials, i.e. on moving away from the neighbourhood of = where 
perturbative calculations are performed for Yang-Mills theories, one arrives at a sufficiently 
large magnitude of for which there exist solutions of Eq. (19) with e — 0, i.e. Eq. 
(18) is satisfied. On further increasing the magnitude of the potentials one finds negative 
eigenvalues, until for an even bigger value of such a magnitude one finds again a solution of 
Eq. (19) with e = 0. 

Note that Eq. (20) is, at this stage, an hyperbolic equation having non-vanishing solutions 
also in the limiting Abelian case; a way to get rid of non- uniqueness of the solutions of the 
equation expressing preservation of the gauge-fixing is to perform a Wick rotation and 
hence consider an Euclidean space. In such a space one can understand the behaviour of 
the solutions of Eq. (20) by exploiting index theory and the concept of spectral flow (see 
Ref. [15] and references therein). 

One can thus imagine that the potentials solving Eq. (18) divide the space of potentials 
into regions, in each of which Eq. (20) has a certain number of eigenvalues, corresponding 
to bound states of the ghost field ipP. In correspondence to the surfaces defined by the 
potentials satisfying Eq. (18), known as Gribov horizons, which separate the various Gribov 
regions, there exist massless ghost states. The first Gribov region C has no ghost bound 
states, the second Gribov region C\ has one ghost bound state, the k-th Gribov region Ck 
has (k — 1) ghost bound states. 

Since the potentials occurring in the Gribov regions C n , n = 1,2,3... are Gribov copies 
of the configurations occurring in Co, as Gribov himself did show in his original paper, the 
functional integration should be restricted to C . Two kinds of Gribov copies can be found, 
i.e. those obtained from the potentials in C from big gauge transformations, resulting 



(20) 
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from topological effects previously described, as well as equivalent copies of the potentials 
in Co that exist within Co itself. While the former can be simply ignored in perturbative 
calculations, the latter make it necessary to use a greater care. To avoid overcounting copies 
in the functional integral, it is necessary to further restrict integration to a domain known 
as the fundamental modular region [16], the boundary of which is studied in Refs. [16-18], 
and within which the Gribov ambiguity no longer occurs [19]. 

In the original analysis by Gribov [14], the restriction of the integration region in config- 
uration space has the effect of removing the infrared singularity of perturbative theory and 
leads to a linear increase of the interaction among non-Abelian charges at large distances, 
i.e. a possible mechanism of quark confinement. 



IV. GRIBOV EQUATION 

Before further discussing the implications of the Gribov analysis, it is appropriate to 
present a simple example, used by Gribov himself, of the occurrence of such a phenomenon. 
An explicit solution of the equation expressing preservation of the Coulomb gauge-fixing 
condition at = under gauge transformations: 

d r ((d r D)D- 1 )=0, (21) 

can be found for the SU (2) group, which can be easily parametrized because the exponential 
map on its Lie algebra is in one-to-one correspondence with the interior of the smallest circle 
where all points map to —I. Thus, every element with the exception of minus the identity 
has a unique representation U : 

U = e iTt ^ = I cos \n\ + in ■ a sin \n\ , (22) 

where n = n/\n\ and n ■ n < ir. To obtain an element U(x) of the gauge group, one has 
simply to make n into a function of x. Let us now assume, to obtain radial symmetry and 
the correct asymptotic behaviour, that [20] 

n=^-(y,x,z), (23) 
r 

where the exchange of x and y is intended, and let us pass to spherical polar coordinates, 
obtaining 

U(r, 9,<P)=[ (24) 
\ — e ^sqSu c u - \ces w J 
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where s a and c a denote sin a and cos a, respectively. To obtain Eq. (24) one has to consider 
the Jacobian matrix for the coordinate transformation and the standard row x column rule 
for the scalar product h • a. On using the familar expressions of the divergence and gradient 
operators in polar coordinates, and the further change of coordinate t = log(r), Eq. (21) on 
U yields [20] 



U e + \uj) C/ t + U t U} + U e Ul + \u^j\ = 0, (25) 

S 9 / S 9 



where subscripts denote derivative with respect to the variable expressed. In the light of 
Eq. (24), Eq. (25) takes the form 



ic e e lv s e 



-e lv s e -ice 



(u + Co — 2 sin oj coscjj = 0, (26) 



where a dot denotes derivative with respect to t. On passing to the variable u; = 2u, and 
leaving aside the 2x2 matrix in Eq. (26) (since it is independent of u), one finds eventually 
a differential equation describing a damped pendulum in a constant gravitational field, with 
a; being the angle from the point of unstable equilibrium. Since U — I at the origin, one 
finds uj = 2kir, where k is an integer. Moreover, since u; is determined up to a multiple of 
27r, there is no loss of generality in taking k — 0. One therefore obtains what is called the 
Gribov pendulum equation 

d 2 d \ 

— + --2sin 57 = 0, (27) 



dt 2 dt 

with the following boundary (or asymptotic) condition: 

lim uj{t) = 0. (28) 

There exist three solutions with the initial condition (28): either the pendulum never falls, 
or it falls anti-clockwise, or instead clockwise. 

An analogous treatment of the Gribov pendulum equation, which however does not rely 
upon the matrix representation, can be found in Ref. [21]. At a deeper level, more than a 
solution is found because the whole of configuration space is being explored. On restricting 
the analysis to a bounded region, as it always occurs in perturbation theory, Gribov copies 
are instead ignored. 
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V. GRIBOV AMBIGUITY AND GENERAL RELATIVITY 



In Sec. II, on discussing the asymptotic conditions for gauge theories, we focused our 
attention on Yang-Mills theories. Actually, the manifestly covariant formalism can treat 
gravitation as well in the same language used for non-Abelian gauge theories. It may be 
interesting to investigate a possible occurrence of the analogue of Gribov ambiguity in general 
relativity. Several ways to express Einstein's theory as a gauge theory have been proposed in 
the literature; the matter is that one can treat this theory just in the connection formalism 
adopted in the introduction for any dynamical system. 

General relativity can be mathematically formulated starting from a pseudo-Riemannian 
4-manifold M. endowed with an atlas of coordinate charts. Symmetries usually considered 
are transformations of the metric tensor induced by general coordinate transformations, 
hence the passive point of view of the diffeomorphism group is preferred to the active one. 
The biggest group Q of passive dynamical symmetries of Einstein's equations is the group 
of transformations reading as [22] 

It is larger than the general coordinate transformation group, usually considered, which is 
a non-normal subgroup of Q [23]. 

In a gauge theory of gravitation one can adopt a principal fibre bundle which has, as 
base space, the union of equivalence classes of all metric tensors, solutions of Einstein's 
equations, generated by passive diffeomorphisms [23] . Despite previous considerations about 
a more general symmetry group involved, such a choice of principal bundle is sufficient to 
characterize the theory. The quotient group where dynamics is reduced, in fact, is the 
same for the three symmetry groups: the passive diffeomorphisms, the active ones or Q. 
The elements of this quotient group are the gauge orbits of general relativity. It is natural 
to investigate topological properties in this case, seeking for the existence of global cross- 
sections, whose lack is found in Yang-Mills theory. 

The differences from the latter theory result from the nature of the diffeomorphism group, 
which is a non-analytic infinite-dimensional Lie group, while Yang-Mills symmetry groups 
are Baker-Campbell-Hausdorff groups [24]. The Gribov ambiguity arises from the non- 
trivial topology of gauge orbits in Yang-Mills theory; in general relativity their topology is 
even less regular, so that the problem becomes more involved as well as more fundamental. 
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Mathematical complications arise from the nature of the diffeomorphism group. Yang- 
Mills theory is a field theory on a background space-time, and the action of the gauge group 
in that case is on a inner space. In general relativity instead, the action of the gauge group 
is an extension to tensors over space-time of the diffeomorphisms of space-time itself [25]. 

It is very difficult, in the manifestly covariant configuration space approach, to face the 
Gribov problem in general relativity. Indeed, one would need tools to make a clean separation 
between gauge variables and a basis of gauge-invariant observables. In other words, it is not 
simple to find suitable coordinates (I A , x a ) as in Sec. I. The Hamiltonian formulation has, 
at least locally, a natural tool for such a separation, i.e. the Shanmugadhasan canonical 
transformations [26]. In this approach, the singularity of Lagrangians, both in particle 
physics and in general relativity, makes it necessary to use the Dirac-Bergmann theory [27] 
of constrained Hamiltonian systems. Since it is the analogue of the space of gauge orbits in 
the global approach, only the constraint submanifold of phase space has physical relevance. 

In phase space there are as many arbitrary Hamiltonian gauge variables as first-class 
constraints, which are the generators of the Hamiltonian gauge transformations. They may 
be used as the I A coordinates of Sec. I, hence they determine a coordinatization of the 
gauge orbits inside the constraint submanifold. To obtain the reduced phase space, one has 
to add as many gauge-fixing constraints as first-class ones. As a consequence of imposing 
gauge-fixing constraints, first-class constraints are turned into the second class and therefore 
treated with the help of Dirac brackets (such constraints are then strongly vanishing). 

The Dirac observables, which in general can be non-local, only give a coordinatization 
of the physically relevant space. On the other hand, it is the analysis of its topological 
properties that makes it possible to say whether a given dynamical system (with constraints) 
admits a subfamily of globally defined Shanmugadhasan canonical transformations. The 
existence of such a subfamily is the same as the one of a global cross-section, from the point 
of view of differential geometry: it means that the system admits preferred global separations 
between gauge and observable degrees of freedom. In other words, in the Hamiltonian 
approach the existence of globally defined Shanmugadhasan canonical transformations makes 
it possible to avoid the Gribov ambiguity. 

Actually, the existence of these transformations has not yet been proved. Indeed, impor- 
tant constraints like the Yang-Mills Gauss law and the ADM supermomentum constraints 
[25] are partial differential equations of ellyptic type, hence they may admit zero modes 
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according to the choice of function space. One has to solve the same kind of problems 
previously treated about the ghost operator in the global approach. 

To obtain, at least as a first approximation, the main non-topological properties of a 
system, the Hamiltonian approach needs to avoid the analogue of Gribov ambiguity in 
general relativity. In fact, in this formalism one assumes that all fields have to belong to 
suitably weighted Sobolev spaces so that the allowed spacelike hypersurfaces are Riemannian 
3-manifolds without Killing vectors. 

We may conclude that in general relativity one has to face the analogue of the Gribov 
problem; in this theory it is more difficult to arrive at some conclusion by virtue of highly 
non-trivial topology of the diffeomorphism group. 

VI. RECENT DEVELOPMENTS IN CONTINUUM GAUGE THEORIES 

We have seen in Sec. II that, to control the asymptotic behaviour of gauge potentials, 
one has to impose appropriate boundary conditions. The choice of these conditions rules the 
occurrence of Gribov copies. However, one could make a distinction between two different 
cases of ambiguity of gauge potentials, according to which kind of boundary conditions are 
imposed. In this way one can have a weak Gribov problem and a strong one [28]. 

When there exist always some configurations which have copies, we can talk of the weak 
Gribov problem. This phenomenon, which is found in any regular gauge, occurs even on a 
compact space and inside a given topological sector of configuration space. 

When there are copies of the vanishing configuration (A^ = 0), hence pure gauge poten- 
tials satisfying the gauge condition, we can talk of the strong Gribov problem. On choosing 
weak-decay boundary conditions, we might find the occurrence of this phenomenon in an 
infinite space. On the other hand, it is not customary inside a given topological sector in 
a compact space; the same occurs on choosing strong decay conditions in a non-compact 
space. 

In terms of functional integration the weak problem is less harmful than the strong one. 
In fact, one might view the strong Gribov ambiguity as a lack of strict positivity for the 
functional measure, and the weak one as a lack of monotonicity. For smooth configurations 
on a smooth compact four-dimensional space, the strong Gribov problem is not expected. 

We have seen that boundary conditions are fundamental in dealing with the Gribov am- 
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biguity. Indeed, appropriate boundary conditions (or compactification) may be imposed to 
avoid the strong Gribov problem, as is usually done in the constructive study of the ultra- 
violet limit of non-Abelian gauge theories to avoid the infrared problem. The weak Gribov 
phenomenon, instead, cannot be avoided [28]. Anyway, in perturbative gauge theories, a 
topological analysis shows that the Gribov problem is irrelevant [29]. 

One could manage to get rid of the Gribov ambiguity by employing different theories, that 
generalize the usual gauge theories. For generalized connections [30], the Gribov problem is 
completely irrelevant for the calculation of functional integrals, hence for the (out | in) ampli- 
tude, on assuming absolute continuity of the considered measure on the space of Ashtekar's 
generalized connections [31] with respect to the induced Haar measure. 

Up to this limitation, the Gribov phenomenon is still present with generalized connec- 
tions, but it is not a problem. Indeed, we have found that topological non-triviality of the 
configuration space is responsible for the occurrence of Gribov ambiguity; this non-triviality 
is concentrated on a zero-measure subset [30] in this formalism. 

Among several techniques proposed to take care of Gribov copies, one might use a quanti- 
zation which does not require gauge fixing [32] such as a generalized Gupta-Bleuler method 
[33], one can modify either the gauge- fixing procedure [34] or the functional integration [35], 
otherwise using non-local variables such that gauge fixing is no longer required [36]. 

Another recently proposed way to get rid of the Gribov ambiguity is a restriction on the 
norm of the (out | in) amplitude, not considering its phase. In this way Yang-Mills theory 
would be free from ambiguities if used for the description of observables [37]. 

The most promising way to avoid the Gribov ambiguity in QCD is, as far as we can see, 
the stochastic quantization approach [38], in which, although there are Gribov copies, they 
have no influence on expectation values. This approach determines an Euclidean probability 
distribution directly in configuration space, i.e. the space of gauge potentials, without 
reduction to the orbit space. 

The first Gribov region may be characterized as the set of relative minima with respect 
to local gauge transformations of the minimizing functional -Fa[-D] = || D ^4|| 2 , analogue of 
the ghost operator JF of Sec. III. The fundamental modular region A, instead, may be 
characterized as the set of absolute minima. The latter is free of Gribov copies, apart from 
the identification of gauge-equivalent points on the boundary <9A, and may be identified with 
the gauge orbit space. 
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One can clearly say that in continuum gauge theories the Gribov ambiguity is still an open 
problem at non-perturbative level. In this respect, one should here mention an important 
result of Neuberger [39], according to which the action used in BRST quantization may 
engender vanishing partition functions, by virtue of Gribov copies which cancel each other's 
effects, and the same may occur for expectation values of physical operators. 

Apart from future developments of the previous methods or other ones to get rid of the 
Gribov phenomenon, for the time being some conclusions may be outlined about implications 
of the existence of Gribov copies in continuum gauge theories. 

First, in the high-energy sector of non-Abelian gauge theory, the restricted domain of 
integration in the functional amplitude has surprising consequences: they contradict the 
expected standard behaviour of the theory [28]. Moreover, the Gribov problem leads to a 
cutoff on frequencies stronger than the expected one, resulting from standard perturbative 
asymptotic freedom [28]. 

In the infrared sector, Gribov himself in his original paper argued that this ambiguity 
leads to an effective infrared cutoff; in fact, he related this behaviour with confinement. 
Last, we may conclude by saying that in non-Abelian theory the Gribov phenomenon seems 
to be an obstruction, at non-perturbative level, to the coexistence of too many energy scales 
[28]. 

VII. THE CASE OF LATTICE GAUGE THEORIES 

In the formulation of lattice gauge theories [40], the gauge-fixing procedure and the 
generation of Gribov copies, obtained via Monte Carlo simulations and, in particular, via 
random gauge rotations [41], are by now well-established techniques, although they have 
been initially ignored since it is not strictly necessary to fix a supplementary (or gauge) 
condition in such theories. The gauge-fixing procedure, however, is particularly convenient 
in some circumstances, and these have led to an investigation of the Gribov problem in this 
framework as well [42]. 

The regularization provided by a lattice [43], in fact, makes the gauge group compact, so 
that the Gibbs average of any gauge-invariant quantity is well-defined. However, asymptotic 
freedom causes the continuum limit to be the weak coupling limit, and such an expansion 
requires gauge fixing. Moreover, to extract non-perturbative results from Monte Carlo sim- 
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ulations, one can use the evaluation of quark/gluon matrix elements, which again requires 
gauge fixing. The same procedure is also used in smearing techniques. 

A brief discussion of the lattice theory formalism may be useful to understand why Gribov 
copies also occur in lattice gauge theory. One can expect their presence by virtue of the 
Killingback analysis [44], which is analogous to the investigation of Singer, but is made for 
Euclidean gauge theories with periodic boundary conditions. 

An action describing fermions and gauge fields on the lattice is obtained by replacing 
differentials with finite differences and constructing a suitable covariant derivative. The 
theory without fermions is known as "pure gauge theory". The pure gauge action on the 
lattice is usually taken to be [45] 

s = PJ2( 1 -w ReTlUp )> (29) 

p 

where N c = 2 for SU(2), and the plaquette variable U p is the ordered product of four link 
variables in a square. The link variables U^ are the degrees of freedom in the theory, and 
are elements of the gauge group. A link between two nearest-neighbour sites i and j of 
the lattice is associated with each link variable U^. Such a variable is also often written as 
Un(x), where x denotes the site of the link and \x the direction. 

Elements of the gauge group are the gauge transformation variables too, which live on 
the lattice sites themselves. A gauge transformation acts on a link variable as follows: 

U^DiUijD' 1 . (30) 

On using the fundamental representation for gauge fields, the action (29) is known as 
the Wilson action. The first term in the expression is a constant and is added in order to 
reduce the action to the Yang-Mills one in the classical continuum limit. This limit is taken 
by letting the lattice spacing go to zero after identifying 

' x J + x 1 



= exp 



ix>-x%A" 



2 

Gauge fixing on the lattice requires to minimize the functional 



(31) 



F L [D] = -J2^^D{x)U^{x)D\x + fi), (32) 



as a function of all gauge transformations D(x). The sum over /i runs from 1 to 3 for the 
Coulomb gauge, and from 1 to 4 for the Landau gauge. Note that this function is the lattice 
analogue of the functional F A [D] = \ \ D A\\ 2 . 
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On expanding U^(x) according to the formal identity [46] 



U„(x) = exp(-a A^x)) = 1 -aA^x) + a 2 Al(x) + ..., 



(33) 



where a is the lattice spacing, one gets 



F L [1] = - Tr (! + a2A K x )) = const - + a 2 F A {l], 



(34) 



where use has been made of the identity Tr A^ = 0. Hence, it is clear that minimizing F L [D] 
is the same as minimizing Fa[D], in the classical continuum limit. 

The functional Fl[D] can be seen as the Hamiltonian of a spin glass system [47]. Such 
a system can behave in a highly chaotic way, hence giving rise to a large number of local 
minima. Indeed, such local minima, which correspond to Gribov copies in the continuum 
gauge theory, appear on the lattice, both for Abelian and non-Abelian fields, in the Coulomb 
gauge as well as in the Landau gauge. The first numerical evidence for lattice Gribov copies 
was presented in the early nineties [41], [47], and [48]. The number of copies in the Coulomb 
gauge was found to be larger than in the Landau gauge [48]. 

The effects on the evaluation of gauge- dependent quantities resulting from the presence of 
Gribov copies are studied by comparing the results for quantities, such as gluon and ghost 
propagators [49], using two different averages: the average considering only the absolute 
minima, which should give the result in the minimal Landau gauge; and the average consid- 
ering only the first gauge-fixed copy generated for each configuration. If Gribov copies were 
not considered, one would obtain the latter average as the result. 

A very interesting analysis is provided by the evaluation of the axial current renormal- 
ization constant Za [41]. The relevance is here twofold: on the one hand, Za can be 
obtained from chiral Ward identities in two distinct ways, either gauge independent (it 
consists in evaluating matrix elements between hadron states) or gauge-dependent (which 
consists in evaluating matrix elements between quark states). This engenders an explicit 
gauge-invariant estimate of Za-, without any Gribov noise (in that the Gribov copies consist 
of fluctuations giving rise to a background noise), which can be directly contrasted with the 
gauge-dependent estimate affected by Gribov copies. On the other hand, there is also the 
advantage of finding Za by solving an algebraic equation of first degree for each time section 
of the lattice, hence avoiding the systematic errors usually emerging from an exponential fit 
of decay signals. 
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The corresponding analysis shows that Gribov copies have visible effects, that can be 
detected, for example, on looking at the slightly different estimates of Za as a function of 
a parameter. The Gribov fluctuations, however, i.e. those induced from the choice of a 
particular Gribov copy, are small and do not prevail on the statistical uncertainty. In other 
words, numerical effects of Gribov copies on a lattice can be divided into two categories: 
measurement distortion and lattice Gribov noise. 

In most investigations of the influence of Gribov copies on lattice quantities, it has been 
found that Gribov noise is of the same order as the numerical accuracy of the simulations, 
and that it scales down as a pure statistical error. In some particular cases Gribov noise 
seems to be quite large. In maximally Abelian gauge it introduces a clear bias on the number 
of monopoles [50] and on the value of the Abelian string tension [51,52]. 

A typical example of distorsion resulting from Gribov copies is the measurement of photon 
propagator in compact U(l) lattice gauge theories. It has been shown numerically [53-56], 
that the photon propagator in the Coulomb phase is strongly affected by Gribov noise and 
that only averages taken on absolute minima of the minimizing functional reproduce the 
theoretical predictions. Thus, one expects that the same might happen for the gluon and 
ghost propagators. Analogue results were found in pure 577(2) lattice gauge theory in the 
minimal Landau gauge [43]. 

In some circumstances, in lattice too there exist copies which do not seem to be related 
by a change in the winding index or from a singular local gauge transformation; they are 
called non-top ological Gribov copies [54]. 

Gribov copies seem to increase in number with increasing lattice volume; experimental 
data, instead, suggest that some Gribov copies are lost as (3 is increased, at fixed volume. 
This is the behaviour that one would expect: a large lattice volume, in fact, is related with 
a large number of gauge transformation variables D{x). Of course, when F L [D] depends 
on more variables, it will have more local minima. Moreover, small (3 means large coupling 
[46] ; in the continuum limit we may expect Gribov copies in such a situation. To obtain the 
continuum limit of the theory we have to let both parameters go to infinity, hence it appears 
highly plausible that the phenomenon of Gribov copies described by lattice theories should 
occur in the continuum limit as well [47]. 

It seems fair enough to say, however, that even for lattice theories no undisputable results 
exist as yet. Investigation of the Gribov problem on a lattice has shown that ambiguities 
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resulting from the gauge-fixing procedure exist for both Abelian and non-Abelian theories 
[57] in the Coulomb and Landau gauges. For example, although on the one hand there is 
evidence that Gribov copies exist in the SU(3) Landau gauge [47], on the other hand the 
measurement of the £77(3) gluon propagator obtained with the help of numerical simulations 
performed in the Landau gauge does not point out the existence of the Gribov ambiguity 
[43]. 

Anyway, every lattice calculation has to face the Gribov problem, for example by employ- 
ing stochastic gauge fixing, as is the case for lattice evaluation of gluon screening masses 
[58]. Furthermore, on looking at the foundations of non-perturbative QCD, if one takes 
locality and BRST symmetry as guiding principles [59], one again finds the unavoidability 
of Gribov copies, which are instead absent at the price of dealing with non-local field theory 
[10]. 
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